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ABSTRACT 



Solar flares and coronal mass ejections (CMEs), the most catastrophic erup- 
tions in our solar system, have been known to affect terrestrial environments and 
infrastructure. However, because their triggering mechanism is still not suffi- 
ciently understood, our capacity to predict the occurrence of solar eruptions and 
to forecast space weather is substantially hindered. Even though various models 
have been proposed to determine the onset of solar eruptions, the types of mag- 
netic structures capable of triggering these eruptions are still unclear. In this 
study, we solved this problem by systematically surveying the nonlinear dynam- 
ics caused by a wide variety of magnetic structures in terms of three-dimensional 
magnetohydrodynamic simulations. As a result, we determined that two differ- 
ent types of small magnetic structures favor the onset of solar eruptions. These 
structures, which should appear near the magnetic polarity inversion line (PIL), 
include magnetic fluxes reversed to the potential component or the nonpoten- 
tial component of major fleld on the PIL. In addition, we analyzed two large 
flares, the X-class flare on December 13, 2006 and the M-class flare on February 
13, 2011, using imaging data provided by the Hinode satellite, and we demon- 
strated that they conform to the simulation predictions. These results suggest 
that forecasting of solar eruptions is possible with sophisticated observation of a 
solar magnetic fleld, although the lead time must be limited by the time scale of 
changes in the small magnetic structures. 

Subject headings: Sun: flares — Sun: coronal mass ejections (CMEs) — Sun: magnetic 
topology — Sun: activity — Magnetic reconnection — Magnetohydrodynamics 
(MHD) 
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Introduction 



Solar flares and coronal mass ejections (CMEs), the most catastrophic eruptions 
in our solar system, are believed to occur through a type of storage-and-release 
phenomena, in which the free ener gy stored as a nonpotential magnetic field in t he 



solar corona is abruptly unleashed ( jPriest fc Forbes 



200d; 



Shibata fc Magara 



20111). It 



is a promising hypothesis that the onset of solar eruptions is related to some typical 
magnetic structures. In fact, many studies have emphasized a relationship between the 



occurrence of so. 



shear ( 1LowII1977 



ar eruptions and various magnetic properties; e.g. strong magnetic 



coronal magnetic field ( IRust fc Kumar 



Hagvard. et al.l ll984l: Mikic fc Linkerlll994l: iKusano et al.lll995l: IZhang 



2001 ) , reversed magnetic shear (IKusano et al.l 12004 : IWang et al.l l2004l: IChandra et al. 



2010l : iRomano fc Zuccarellol l201ll: IVemareddv et al.ll2012r). sigmoidal structure of t he 



cancellation on t 



Green et al 



Amari et al. 




le photosphere (Ivan Ballegooiien fc Martens 



Canfield et a1.lh999l: iMoore et allbomh. fl ux 



198f 
11 



Zhang et all 1200 



IT' 



2003a 



co nverging foot point motion (Jnhester et al...l992i : iForbes fc Priestlll995 



Cheng et al. 



emerging magnetic fiuxes (jHevvaerts. Priest fc Rust 



Fevnman fc Martin 



Walla 



,ce et al. 



2OIOI: 



995 



2010h. the sharp gradient of magnetic field (ISchriiver 



Chen fc Shiba 



Moore et al. 



topol ogies (lAntiochos et al 



rope ( IForbes fc Priest 



1995 



1 



2011 



2OOOI: 



1977: 



Moore fc Roumeliotis 



Fan fc Gi 



Archontis fc Hood 



Sterhng fc Moor 



Torok fc Kliem 



2005 



2iM 



Chen 



3Sonl 12004 
2OI2L multip olar 



Chifor. et al. 



Willi 



19961: 



ams et al 



narro w magnetic lanes betwee n major sunspots f^Zirin fc Wang 



2007h. topological complexitv fISchmieder et al. 



( lAbramenko fc Yurchyshyn 



'an fc 



199^: 



Gibson 



20051), fiu x 



20071), 



Kubo et al. 



I994I ) , interm ittency and mu ltifractality 



2010), and double loop structure ( IHanaokal 119971 ). 



However, although some observations suggested that these properties tend to correlate 
with the solar eruption productivity, the triggering mechanism for solar eruptions is still 
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debated ( IForbes 



solar eruptions ( ILeka fc Barnes 



2000r). and we do no t 



2007 



yet well understan d what determines the onset of 



Welsch et al. 



20111 ). It substantially hinders our 



capacity to predict the occurrence of solar eruptions and to forecast space weather. The 
objective of this paper is to clarify the types of magnetic field structures that are capable 
of generating solar eruptions to gain a better understanding of their triggering mechanism. 

Because the magnetic structures of eruption-productive active regions (ARs) are highly 
complex, it is difficult to determine the key elements of the magnetic field present at the 
onset of eruptions from the imaging data. In this study, we compensate for this difficulty 
by systematically surveying the nonlinear dynamics caused by a wide variety of magnetic 
structures, which consists of a large scale force- free field and a small scale bipole field, 
in terms of three-dimensional magnetohydrodynamics (MHD) simulations. As a result, 
we show that two typical magnetic structures favor the onset of solar eruptions, and we 
demonstrate that th e simulation results are in good agreement with the observations with 



the Hinode satellite ( IKosugi et al. 



mi). 



The paper is organized in the following manner. The model and results of simulations 
are described in Sections 2 and 3, respectively. On the basis of the simulation study, we 
propose the two scenarios for the triggering process of solar eruptions. In Section 4, we 
show that two large fiares observed by Hinode, the X-class fiare on December 13, 2006 and 
the M-class fiare on February 13, 2011, conform to the predictions determined through 
simulation. Finally, in Section 5, we discuss the topological features of the magnetic field in 
the triggering phase of solar eruptions as well as the predictability of these eruptions on the 
basis of this paper's conclusion. 
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2. Simulation Model 

The simulation box includes a rectangular domain (— 3L, — 0.75L, 0) < {x,y,z) < 
(3L, 0.75L, 3L), which corresponds to the solar atmosphere above the photosphere within 
an AR. L is the typical length and the coordinate z represents the altitude from the 
photosphere. The x axis [y — z — Qi) corresponds to the magnetic polarity inversion line 
(PIL), on which the vertical component of magnetic field Bz changes sign. 

To model the prefiare state, the initial magnetic field is given by the linear force-free 

field, 

^ ak~^ cos ky ^ 
^init{oi] z) = So -Kk~^ cos ky exp(-i^^z), 
y sin ky j 

where k = 27r/L, K ^ [k? — a^Y^^^, and is a constant. This field is characterized by the 
scalar parameter a and shear angle 9^ = tan~^(— S^^/S^) = tan~^ a/ which is defined as 
the azimuthal rotation of the magnetic field with respect to the potential magnetic field 
Bimt(0; 0, 0) = (0, -So, 0), as shown in Fig. 1. 

After the simulation starts, we quickly inject a small bipole magnetic field into the 
force-free field ^init from the bottom boundary in order to form the partially disturbed 
magnetic structure illustrated in Fig. 1. This process is performed by prescribing the 
boundary condition of electric field, which is generated by the constant ascending motion 
of a synthetic magnetic torus located virtually below the simulation box. 

The ascending torus forms a sphere of radius rg, within which there is only the toroidal 
magnetic field Bg of constant intensity. The major axis of the torus is on the x — y plane 
and its center is initially at a point (0,i/e, — ^g). The torus ascends with constant velocity 
Ve for only the period < t < re(< r^jv^. The electric field Bg x Vg and velocity Vg are 
imposed within the cross-section of the ascending torus on the bottom plane until t = Tg. 
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The injected small scale field is characterized by the azimuthal orientation Lp^ of Bg at the 
top of torus, the offset of torus center y^^ and the total amount of injected fiux that is a 
function of Sg, Tg, and v^Te (see Fig. 1). 

Although this m odel corresponds to the kinematic model of emerging fiux 



(IFan fc Gibson 



20041 ). we adopt it as a way to dynamically form a variety of magnetic 
structures rather than as the model of fiux emerging. In fact, because the speed of small 
magnetic field injection is about two orders of magnitude faster than the realistic speed 
of emerging fiux (the order of 1 km/s), the physically meaningful process of simulation is 
after the fiux injection terminates. Therefore, we focus mainly on the dynamics caused by 
magnetic structures formed as the result of fiux injection rather than the dynamics of fiux 
emerging. Which process in fiux emerging or horizontal motion on the solar surface is more 
eflScient for triggering the solar eruptions is beyond the scope of this paper, although we 
partially discuss this issue in Section 3. 

The length, magnetic field, and time are described by non-dimensional units, L, 5o, 
and the Alfven transit time {ta = L/Va)^ respectively, where Va is the Alfven speed 
Bq/ -^/LionrUi for ion mass m^, plasma number density n and vacuum permeability /xq. In a 
typical AR {Bo - 0.05 T, n - 10^^ m-^ L - 20 Mm), ta ^ 0.4 s and Va ^ 50 Mm/s. 

We adopted the zero-beta model for the simulation, in which plasma pressure and 
variations in plasma density are assumed to be less effective. This model is applicable when 
the plasma beta (i.e., the ratio of thermal to magnetic pressure) is significantly less than 
unity. The coronal plasma within an AR in the prefiare phase, which is the focus of this 
paper, is thus covered b y this model. The governing equ ations are the same as those used 



in our previous study ( IKusano et al. 



2004: 



Kusano 



20051 ). The electrical resistivity r] and 



viscosity u are initially constant (10 ^ and 5 x 10 ^ in non-dimensional units, respectively). 
However, rj increases if the electric current density J exceeds a critical Jg, as adopted from 
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Eq. 9 of iKusano et al.l ( 120041 ). where the enhanced resistivity ryi = 5 x 10 and critical 



electric density Jc = 50. 

The spatial differential operator is approximated by the second-order accurate finite 
difference method using three-grid-point stencils, and the temporal integration was 
performed using the Runge-Kutta-Gill method with fourth-order accuracy. The grid 
number included in the simulation box is 1024 x 256 x 512 for each dimension x, and 

and they are packed near the PIL so that the finest grid sizes are Ax = 5.9 x 10~^, 
Ay = 2.9 X 10"^, and Az = 5.9 x 10"^, respectively. 



3. Simulation Results 

We chose the shear angle 6q of large scale force-free field and the azimuthal orientation 
of the small scale injected field as the parameters, and surveyed the parameter space 
(0, 0) < (^0, ^e) < (85°, 360°) with 106 runs. The cases = 0°, 0° < (^e < 180°, = 180°, 
and 180° < (fe < 360° correspond to the small scale field of the right polarity, normal 
shear, opposite polarity (OP), and reversed shear (RS), respectively, compared with the 
preexisting large scale field. 

The simulation results are summarized in Fig. 2, which clearly shows that the kinetic 
energy produced by eruption strongly depends on Oq^ hence a large-scale eruption is possible 
only in strongly shear ed cases (^o > 75°). This result agrees well with the observations 



( IHagyard. et al. 



19841 ). and is logical because a more strongly sheared field stores a greater 
amount of free energy as the nonpotential component of the magnetic field. However, the 
most notable feature of this figure is that the occurrence of solar eruption is sensitive to 
(fie- In particular, the eruption-producing cases, represented by diamonds in Fig. 2, exist 
mainly for Lp^ = 123° to 270°. This result indicates that strong shear alone is not a sufficient 
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condition for the onset of eruptions and that the occurrence of eruptions is governed 
by small magnetic structures. The simulation results predict that the OP and RS-type 
magnetic structures are capable of triggering solar eruptions. 

We detected a clear difference in the morphologies of magnetic fields between cases 
triggered by OP and RS-type configurations, represented in Fig. 2 by pink and blue 
diamonds, respectively. The typical dynamics of eruption caused by the OP-type field are 
explained in the following steps, as shown in Fig. 3: (1) Two sheared magnetic field lines 
rooted near the small bipole of OP (blue tubes in Fig. 3a) are reconnected via the bipole of 
OP to form twisted fiux ropes (green tubes in Figs. 3b-3d), (2) after twisted fiux ropes grow 
(Fig. 3e), they suddenly erupt upward (Fig. 3f), (3) in the strongly sheared case (^o > 70°), 
the eruption of twisted fiux ropes vertically stretches the overlying field lines, and new 
vertical current sheets are generated below them (red surface in Fig. 3f), and (4) magnetic 
reconnection begins on this vertical current sheet and forms a cusp-shaped post-fiare arcade 
along with more twisted ropes, which accelerate the eruption (Fig. 3g). In particular, for 
the most strongly sheared case (^o > 75°), steps (3) and (4) reinforce each other and the 
reconnection region propagates along the PIL (Fig. 3h) until reaching the boundary of the 
simulation box, which corresponds to the outer border of the AR. 

As a result, in a more sheared arcade, greater kinetic energy is produced by the ascent 
of longer fiux ropes. This process is essentially the same as the model of tether cutting with 



Moore fc R.oumeliotisl (119921) . The rapid ascent of the twisted 



emerging fiux reported by 
rope can be attributed to the loss of equilibrium and id eal MHD stability (iForbes fc Priest 



199^; 



Kliem fc Torok 



200^; 



Demoulin fc Aulanier 



mm) 



On the other hand, the RS-type fiux triggers eruptions in a different manner. As seen 
in Fig. 4, the injected bipole field of reversed shear contacts the preexisting sheared fields 
(blue tubes in Figs. 4a-4c) and forms a current sheet (red surface denoted by c in Fig. 4b) 
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on the border between them. Magnetic reconnection slowly proceeds on this current sheet, 
and the sheared field is removed from the center to the sides of RS field region (d and d' in 
Fig. 4b). Because of the reduction of the sheared field in the center, the magnetic arcade 
above the injected small bipole partially collapses into the center, and a vertical current 
sheet (v in Fig. 4d) is generated. As a result, sheared magnetic fields (sky-blue tubes in Fig. 
4d) are reconnected to form twisted fiux ropes rooted on F and F' in Fig. 4e. These twisted 
flux ropes grow (Fig. 4f) and erupt upward (Fig. 4g). Finally, the processes continue in the 
same manner as those in steps (3) and (4) for the OP-type case (Fig. 4h), and the eruption 
develops. 

This triRg;er s cenar io is essentially the same as the reversed-shear flare model proposed 



by iKusa no et al. 



However, it should be noted that they studied the process 
that horizontal flow along the PIL, rather than flux emerging, reverses magnetic shear. 
Therefore, we can conclude that emerging flux is not necessarily needed for the onset of 
eruption, if horizontal motion drives the formation of the magnetic structures favorable to 
the onset of eruptions. 

It should be noted that the OP-type conflguration may cause an eruption of flux 
rope before the onset of flare reconnection, whereas the RS-type conflguration triggers 
flare reconnection before the eruption of flux rope. We refer to these processes as 
"eruption-induced reconnection" and "reconnection-induced eruption," respectively. Figure 
2 indicates that while the eruption-induced reconnection is possible for any value of sheared 
angle Oq^ reconnection-induced eruption can occur only for relatively large ^o- We can easily 
explain this result because some amount of magnetic shear has to pre-exist for reconnection 
of the RS fleld to cause partial collapse of the magnetic arcade. 

In both processes triggered by the OP and RS-type conflgurations, the morphology 
of the magnetic fleld in the latter phase (after the onset of flare reconnection) is common 
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and consistent with the standard CSHKP model of a tw o-ribbon flare (jCarmichael 



Sturrock 



19661 : 



Hirayama 



1974: 



Kopp fc Pneuman 



1964 : 



19761 ). However, the initiation procedure 



differs between the two configurations. It indicates that although both the eruption of fiux 
rope, which may create CMEs, and the reconnection of flares can induce each other, the 
causality between them may be governed by small magnetic structures that trigger the 
eruptions. 



4. Comparison with Observations 



On the basis of our simulations, we determined that the magnetic structures favorable 
to the onset of large solar eruptions consist of the OP or RS-type small magnetic flelds 
and the strongly sheared fleld, as illustrated in Fig. 5. For either type of trigger processes, 
major flares should be preceded by minor reconnection between the preexisting sheared 
fleld and the flux element with a different orientation. These results are consistent with 



the observations of internal reconnection ( IMoore et al. 



200ll ). Therefore, we can anticipate 



preflare brightening related to the precedent reconnection in both OP or RS-type regions as 
a precursor to major eruptions. In fact , we have determined that two major flares observed 



by the Solar Optical Telescope (SOT) f Tsuneta et al. 
prediction. 



20081 ) aboard Hinode conform to this 



4.1. Case of 2006 December 13 Flare 



The flrst event is the X3.4-class flare observed in AR NOAA 10930 at 02:14 UT on 



December 13, 2006 ( IKubo et al. 



20071 ). According to the vector magnetogram observed by 
the SOT/Stokes Spectro Polarimeter (SP) installed on Hinode, the azimuthal magnetic fleld 
in this region was mainly twisted in a clockwise direction by more than 75°. We analyzed 
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the spatiotemporal correlation between the hne-of-sight component of the magnetic field 
and the Ca II H line emission, which were observed by Narrowband Filter Imager (NFI) 
and Broadband Filter Imager (BFI) of SOT, respectively. As evident in Fig. 6a, a small 
iso lated positive pole (denoted by I) slowly grew since 23:00 UT on December 12. According 



to 



KuboetaL 



( 120071 ). this small magnetic flux emerged around UT on December 12 at 
the west side of the positive umbra, and it was elongated by the rotating flows around 
the southern positive umbra. As a consequence, an OP-type conflguration was formed in 
the region denoted by a yellow circle (Fig. 6b), in which the orientation of a small bipole 
(positive in the north and negative in the south) was opposite to the orientation of positive 
and negative major spots, as illustrated in the subset at the lower right corner of Fig. 6b. 
A brightening of the Ca II H line began on the PIL of this OP-type bi-pole 30 min prior to 
the flare onset (Fig. 6b) and extended to the flare ribbons (Figs. 6c-6i). 

In the very early stage of the flare (Figs. 6c and 6e), the bright kernels consisted of 
sheared ribbons (F on the positive pole; F' on the negative pole) and small segments (Pi, 
P2, and P3). Because the chirality of shear ed ribbons is con sistent with the observation 



of negative magnetic helicity in this region ( jPark et al. 



20101 ). the sheared ribbons can be 



interpreted as the feet of twisted flux rope formed by reconnection in the OP-type trigger 
region as seen in Figs. 3b and 3c. Here it should be noted that the feet F and F' in Fig. 
3b corresponds to a mirror image of the observed F and F' in Fig. 6c because the sign of 
preloaded magnetic helicity differs. The sign is negative in this AR and positive in the 
simulation. 

The transient bright segments (Pi, P2, and P3) are consistent with the simulation 
results (Fig. 3), assuming that segments appeared on the intensive electric current sheets 
formed on the top, bottom and lateral side of the twisted rope indicated as Si, S2, and 
S3, respectively, in Figs. 3b-3d. The current sheet Si is formed between the twisted rope 
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and overlying field represented by blue tubes in Fig. 3d. Therefore, if a chromospheric 
emission Pi corresponds with the curren t sheet Si , it imp lies that the lowest part of the 



twisted rope, which forms a bald patch ( ITitov et al. 



I993I ). existed in the chromosphere. 



The fact that the bright segment Pi quickly disappeared at 02:20 UT (Fig. 6e) can be 
attributed to the ascension of the fiux rope beyond the chrom osphere. In the soft X-ray 



images observed by X-Ray Telescope (XRT) (.Golub et al, 



20071 ) installed on Hinode, a faint 



ejection is visible from 02:18:18 to 02:22:18 UT, as indicated by arrows in Figs. 6f and 6h. 
This faint ejection exhibited a bar- like structure along the magnetic neutral line and the 
ejection speed changed from 40 to 100 km/s. It is likely that the ejection observed by XRT 
corresponds to the front of the fiux rope launched from the chromosphere to the corona. 

While the P2 segment was weakened (Figs. 6c, 6e, and 6g), the P3 segment grew 
and was combined with ribbon F (Fig. 6g). In addition, these short-term evolutions of 
the bright segments agree well with the simulation result in which the current sheet S2 is 
weakened; however, S3 grows to the major current sheet below the erupting fiux rope (Figs. 
3d-3f). Finally, the ribbon on the negative pole (F' in Fig. 6e) was extended eastward to 
form ribbon R2 and the ribbon on the positive pole F thickened, as denoted by Ri in Figs. 
6i and 6j. The structure and dynamics of these main two-ribbons are consistent with the 
evolution of the feet of the post-fiare arcade observed in the simulation results (Ri on the 
positive pole and the counterpart of R2' on the negative pole in Figs. 3g and 3h). From this 
excellent agreement with the observation and simulation, we can conclude that the onset of 
this event can be explained by the OP-type trigger model. 



4.2. Case of 2011 February 13 Flare 



The second event, which can be explained by the RS-type trigger model, is the 
M6.6-class fiare observed in AR NOAA 11158 at 17:28 UT on February 13, 2011 (Fig. 7). 



-13- 



The vector magnetogram (SOT/SP) indicates that the magnetic field in this region was also 
extensively twisted (more than 75°) in a counterclockwise direction. This counterclockwise 
twist, which corresponds to the preloaded positive magnetic helicity, is consistent with the 
sheared structure of the fiare ribbons (Fig. 7d). 

Here we emphasize that an inverse S-shaped structure of the PIL, in which positive and 
negative polarity regions edge into each other, was formed in the region denoted by yellow 
circles in Figs. 7a- 7c between 13:15 UT and 15:00 UT. On the center of this region, the PIL 
is running north-south, whereas the overall PIL generally is running east- west. Therefore, 
it is highly likely that magnetic field crossing the inverse S-shaped PIL is reversely directed 
compared with the large scale sheared magnetic field, and that the RS-type configuration is 
formed as illustrated in the subset at the lower right corner of Fig. 7b. 

In fact, the brightening of the Ca II H line (Pi in Fig. 7b) was observed since 15:00 
UT in this area. This prefiare brightening can be attributed to heating of the current sheet 
between the pre-existing sheared field and magnetic field of RS, which corresponds to c in 
Fig. 4b. Immediately after the onset of the fiare, sheared ribbons F and F' appeared as a 
center on the inverse S-shaped PIL (Fig. 7d). The ribbon structure is highly consistent 
with the simulation result such that the feet of the twisted ropes (F and F' in Fig. 4) and 
post-fiare arcades (R and R') combined to create elongated sheared ribbons on both sides of 
the RS field region (Fig. 4h). 

Furthermore, it is worthwhile to note that the prefiare brightening (Pi in Fig. 7) was 
weakened before the onset of the fiare (see Fig. 7c). In the simulation, the electric current 
intensity on the current sheet c was weakened just prior to eruption of twisted ropes (see 
Fig. 4d). This weakening resulted from the reconnection on the current sheet v, because 
the reconnection ejected downward the shear-free field, which prevented contact between 
the pre-existing sheared field and the field of RS. The agreement of the observation and 
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simulation results indicates that the weakening of preflare brightening on the RS-hke region 
could be regarded as a precursor of flare onset. 

From the magnetic structure and evolution of preflare and flare ribbons, we can 
conclude that this event is consistent with the RS-type trigger scenario. The formation 
process of the inverse S-shaped PIL in this region is an interesting subject and we are 
analyzing the cause of the strucutre. The detail results of analysis will be reported 
elsewhere. 



5. Discussion and Conclusion 

The results shown in the previous section indicate that events indeed occurred, which 
are described by the following two scenarios: Scenario 1 involves an eruption-induced 
reconnection process triggered by an OP-type magnetic fleld, and Scenario 2 is a 
reconnection-induced eruption process triggered by an RS-type magnetic fleld. Here, we 
examine the reason for existence of both scenarios. It is explicable in terms of the topology 
of the magnetic fleld. In both scenarios, flare reconnection on a vertical current sheet is 
caused by the diverging flows that remove magnetic flux and plasma from the reconnection 
site. This mechanism corr esponds to the "pull" mode in reconnection experiments 



performed in laboratories ( lYamada 



19991 ). As shown in Fig. 8, two different types of 
divergent flows are able to form thin current sheet v between two sheared magnetic loops 
Bi and B2. One divergent flow is a vertical (upward) flow above the PIL (u in Fig. 8), and 
the other is a horizontal flow along the PIL (h and h'). They correspond to the preflare 
flows that trigger flare reconnection in Scenarios 1 and 2, respectively. Because a divergent 
flow can flow either in the vertical or horizontal direction, two different scenarios exist. 

On the other hand, the converging flow (c and c') can also form a vertical current sheet 



- 15 - 



and may driv e flare reconnec t ion, which corresponds to the "push" mode in reconnection 



experiments (lYamada 



19991 ). In fact, the converging flow on the solar surface is behey ed 



to play the role of a tri^^er in the flu x cancellation (FC) model ( jinhester et al. 



Priest et al. 


1994; 


Amari et al. 


2003a) 



1992 



topology, we can conclude that these three scenarios (Scenarios 1 and 2 proposed in this 
paper as well as the reconnection scenario caused by the converging flow) can all apply to 
the triggering of solar flares. 

The classiflcation of the PIL morphology in flaring sites is a promising technique 
used to examine the applicability of each scenario because each scenario requires different 
structures of PIL as illustrated in Fig. 5. Scenario 1 requires an OP-type disturbance to 
the PIL (Fig. 5a). On the other hand, the inverse S-shaped PIL (RS+ in Fig. 5b) and 
the S-shaped PIL (RS— in Fig. 5b) are the notable signature of the RS structures when 
the large-scale sheared fleld is positive (counterclockwise rotated) and negative (clockwise 



rotated), respective. 



(IKusano et al. 



y, although an RS structure without a disturbed PIL is possible 



20041 ). In fact, such types of the PIL morphology have be en observed in 



severa l flare-productive active regio ns; e.g. ARs NQAA 66 59 (Fig. la of 



(119931 )) and NO A A 9026 (F ig. 4q of 



AR NQAA 9026 (Fig. 4k of 



Kurokawa et al. 



Kurokawa et al 



Zirin fc Wang 



(120021 )) for an OP-type morphology. 



(120021 )) for an RS-type morphology. 



On the other hand, the FC model caused by the converging flow does not require 
any disturbances to the PIL if positive and negative magnetic fluxes cancel each other. 
However, because the emergence of a magnetic element in the flaring site is often associated 



with the cancellation of magnetic flux ( IWang fc Shi 



19931 ) , we should carefully recognize the 



structure of flare triggering. Our analyses suggest that a high cadence (less than a couple of 
minutes) of preflare brightening observation in the chromospheric lines such as Ca II H is a 
powerful tool to recognize preflare processes. The statistical study of various solar eruptions 
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will explain the applicability of each scenario to the onset of solar flare and CMEs. 

Our study suggests that solar flares and CMEs are triggered by magnetic disturbances 
in regions where magnetic shear has stored free energy rather than originating in phenomena 
that inevitably occur as a consequence of free energy storage; i.e., both magnetic shear (or 
magnetic helicity) and "proper" disturbance of the magnetic fleld are necessary for the 
occurrence of solar flares and CMEs. This conclusion implies that the lead time for flare 
and CME forecasting is limited by the time scale of changes in small magnetic elements 
that work as a trigger. Our analysis of the two events suggests a lead time shorter than 
several hours. Longer forecasting can be difficult and may be possible only in a probabilistic 
manner. 

While the simulations indicate that the orientation of a small magnetic element that 
perturbs the sheared magnetic structure is crucial for triggering flares and CMEs, it is 
highly likely that the size and location (displacement from PIL) of magnetic perturbation 
are also important. The results of more advanced simulations to determine the significance 
of these parameters will be published elsewhere. 
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Fig. 1. — Illustration of simulation setup. Curved solid lines with arrows correspond to 
magnetic field lines in the preexisting force-free field and the small scale injected field, which 
rotate ^iid (/:?e, respectively, with respect to the large scale potential field. White and gray 
areas on the bottom surface indicate positive and negative magnetic polarity regions. 
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Fig. 2. — Summary of simulations for = 2, rg = 0.13, Tg = 20, = 6.7 x 10"^, ?/e = on 
parameter space of (pe and ^o- Different marks (squares and diamonds) represent the types of 
dynamics, and contours show the maximum total kinetic energy produced by eruption {E^). 
Squares indicate that no eruption has occurred at the corresponding parameter; diamonds 
indicate the appearance of eruptions at each parameter. Pink and blue diamonds indicate 
eruption-induced reconnection and reconnection-induced eruption processes, respectively. 
The yellow diamond corresponds to a special case in which the potential field collapses 
because of reconnection with the small scale injected field, which exhibits a completely 
antiparallel polarity compared with the initial potential field. Right hand and top subsets 
illustrate the initial sheared field and orientation of injected small bipole field, respectively, 
in which white and black areas indicate positive and negative polarity and arrows represent 
the horizontal component of the magnetic field. 
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Fig. 3. — Simulation result for ip^ = 180° and = 77.5°, in which the OP-type of magnetic 
structure causes the eruption-induced reconnection dynamics. Each subset represent a birds 
eye view (a, c, e-h), top view (b), and enlarged side view (d) of the magnetic field at 
different times. Green tubes represent magnetic field lines with connectivity that differs 
from the initial state. Selected magnetic fields in the initial state and those retaining the 
initial connectivity are plotted by blue tubes in (a) and (d), respectively. Red contours 
correspond to intensive current layers in which |V x B| > 40. Gray scales (white, positive; 
black, negative) on the bottom plane indicate the distribution of the z component of magnetic 
field Bz. Other parameters are = 2, rg = 0.13, Tg = 20, Ve = 6.7 x 10~^, and = 0. 
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Fig. 4. — Simulation result for Lp^ = 270° and 6q = 77.5°, in which the RS-type of magnetic 
structure causes the reconnection-induced eruption dynamics. Each subset represents a 
birds eye view (a, c, and e-h) and enlarged side views (b and d) of the magnetic field at 
different times. Green tubes represent magnetic field lines with connectivity that differs 
from the initial state. Selected magnetic fields in the initial state and those retaining the 
initial connectivity are plotted by blue tubes in (a) and blue and sky-blue tubes in (b-d), 
respectively. Red arrows in (b-d) represent the typical flow directions. The format and other 
parameters are same as those described in Fig. 3. 
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Fig. 5. — Schematic diagram of top view of (a) OP and (b) RS-type magnetic structures. 
Arrows indicate the large scale sheared magnetic field, and black and white parts correspond 
to the negative and positive polarity regions. The chirality of PIL disturbance in the RS- 
type structure is changed depending on the sign of magnetic shear; RS+ for positive shear 
(counterclockwise rotation) and RS— for negative shear (clockwise rotation). 
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Fig. 6. — Onset processes for solar flares in AR NOAA 10930 observed by Hinode/SOT 
(a-c, e, g, and i), and soft X-ray images obtained with the Hinode/XRT (d, f, h, and j). 
Gray scale in the SOT data represents the Stokes- V/I signal of the Fe I 6302 A line, which 
indicates circular polarization corresponding to the line-of-sight component of the magnetic 
fleld (black, negative; white, positive). Red contours indicate the emission of the Ca II H 
line. Green lines on the SOT data and gray lines on the XRT data represent the polarity 
inversion line (PIL) in the central region. Observation times are shown for each subset. 
Subset at the lower right corner of (b) illustrates the typical OP-type magnetic structure 
with clockwise rotating sheared field. 
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Fig. 7. — Onset processes for solar flares in AR NOAA 11158 observed by Hinode/SOT. Gray 
scale (black, negative; white, positive) represents the Stokes- V/I signal in the Na I 5896 A 
line. Red contours indicate the emission of the Ca II H line, and green lines represent the 
polarity inversion line (PIL). Observation times are shown for each subset. Subset at the 
lower right corner of (b) illustrates the typical RS structure with counterclockwise rotating 
sheared fleld. 
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Fig. 8. — Illustration of the three different flows u, h-h', and c-c', which may cause the 
formation of a vertical current sheet v and magnetic reconnection between two sheared 
magnetic loops Bi and B2 on the polarity inversion line (PIL). 



